Non-exclusion statistics: a generalization of Bose-Einstein's principle 
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By constructing the super-particle representation of the free boson gas, we propose a new statistics 
in which the particles are non-exclusive. This statistics can be considered as a generalization of 
Bose-Einstein's. The possible condensation of this statistical system is studied. It is found that the 
chemical potential below the condensation temperature is linearly proportional to the temperature 
rather than a constant. With an proper choice of the exclusion factors 7;, Hadane-Wu's fractional 
statistics is retrieved in this representation. 



Bose-Einstein statistics and Fermi-Dirac statistics are 
two paradigms of quantum statistical mechanics. The 
former is non-exclusive, i.e., an arbitrary number of parti- 
cles (bosons) are allowed to occupy a single-particle quan- 
tum state, while the latter is exclusive and in maximum 
only one particle (fermion) is allowed to occupy a single- 
particle quantum state. In a condensed matter, however, 
the interaction may change the statistics of the quasi- 
particles drastically. A well known example is the one- 
dimensional hard-core boson lattice, which is identical to 
a free fermion lattice via the Jordan- Wigner transforma- 
tion. That means the infinite on-site repulsion changes 
the statistics of the bosons to a Fermi-Dirac one. In 
fact, any small repulsive interaction in a one-dimensional 
boson system may drive it to behave as a Fermi sys- 
tem (in the sense that the system has a Fermi surface), 
as shown by Yang and Yang in an integrable model. 
Motivated by the study of spinons in one-dimensional 
antiferromagnets, J&Q] Haldane formulated a general ex- 
clusion statistics S without specific reference to spatial 
dimension. With Haldane's definition, the maximum oc- 
cupation number of particles in a single-particle quantum 
state is just the inversion of the statistical factor g, as 
shown by Wu. f| Therefore this definition can be con- 
sidered as a generalization of Pauli's exclusion principle. 
When applied to fractional quantum Hall effect, it coin- 
cides with the standard 2D anyon definition JtJ in terms 
of the braiding of particle trajectories, though in gen- 
eral they are different. There are indeed some systems in 
which the interactions among the bare particles can be 
completely absorbed into the quasi-particles with which 
the systems behave as an ideal statistical ones. Concrete 
examples are the Calogero-Sutherland model ^,|) and 
the Haldane-Shastry spin chain. ||^] 

Up to date, a successful Bose liquid theory correspond- 
ing to Landau's Fermi liquid theory has not been estab- 
lished. The main difficulty lies in the absence of an en- 
ergy cutoff (counterpart of the Fermi energy in a fermion 
system). Though an arbitrarily small repulsive interac- 
tion may drive a one-dimensional boson system to be- 
have as a non-ideal exclusion system or a Lut- 
tinger liquid, pl| the superfluidity of the liquid 4 He and 



the condensation of some alklai atoms strongly indicate 
that the interactions of bosons may not induce an ex- 
tremely exclusion in three dimensions. In this letter, we 
formulate an ideal non-exclusion statistics beyond that of 
Bose-Einstein. Comparing to Haldane's exclusion statis- 
tics, this new statistics can be considered as a general- 
ization of Bose-Einstein's principle. The deviation from 
Bose-Einstein statistics can be due to some interaction 
which is supposed to be absorbed into the quasi-particles 
completely. 

For a free boson system, an arbitrary number of parti- 
cles occupying a single-particle quantum state is allowed. 
In addition, the physical quantities of these particles are 
additive, i.e., if n particles are in the quantum state la- 
beled by the momentum k with single-particle energy 
e(k), their total momentum and energy take the value 
of nk and ne(fc), respectively. Therefore, I bosons with 
mass M in a single-particle state k can be treated as 
a super-particle of mass IM, momentum Ik and energy 
le(k). We call it simply a rank-/ super-particle. Ob- 
viously, a rank-/ super-particle can also be treated as 
the combination of some rank-/ Q super-particles with 
J2 a la = I- If we treat all types of super-particles in- 
dependently, the many-particle quantum states must be 
over counted. Hence an intrinsic exclusion among the 
super-particles defined in the same single-particle quan- 
tum state should be taken into account to give a cor- 
rect statistics. In the thermodynamic limit, we define 
the density of rank-/ super-particles in the single-particle 
phase space as pi(k) and the corresponding hole density 
as pHk), respectively. The density of particles (rank-1) 
is given by n(k) — ^2iZilpi(k). The intrinsic exclusion 
of the super-particles is expressed formally as 



flmpm(k), 



(1) 



where p® denotes the density of rank-/ holes in the ab- 
sence of particles and fi m are the intrinsic exclusion fac- 
tors. Since the volume of a boson in the phase space 
does not depend on its mass and we measure the den- 
sities in the phase space of rank-1 particles, we readily 
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have p° = Ipl with p\ = l/(2irfr) d in d spatial dimensions 
as usual. The linear dependence of p 1 ^ and p m is up to 
now still an assumption. We note that if Eq.(l) gives a 
correct description of a free boson system, the exclusion 
factors fim ought to be dimension independent. To show 
them explicitly, let us consider a one-dimensional system 
consisting of N interacting bosons with (^-potential. The 
Hamiltonian reads 



A' 



H : 



T — 

3=1 3 



2c S(xi — xj ). 



(2) 



l<3 



This model was shown to be exactly solvable [fL0[ for ar- 
bitrary real c and for c — ► 0, it is reduced to a free boson 
system. The eigenstates of Eq.(2) are described by N mo- 
menta {kj} which are determined by the so-called Bcthc 
ansatz equations 1 1 1 



/ = ]T / *[ e $ - APi$)dk-T / s(k)dk, (7) 
i=i J J 

where p, denotes the chemical potential, e(k) is the energy 
of a single particle and s(k) is the entropy density which 
reads |Q 

oo 



1=1 



- Pl (k) In Pl (k) - p?(k) In p> l (k)}. (8) 



By minimizing the thermal potential via Spi(k) with 
Eq.(l) and Eq.(6), we obtain 



(9) 



N 

■m— 1 J 



(3) 



with the eigenenergy of Eq.(2) as E = kj, where L is the 
length of the system. For c — > CP, it has been demon- 
strated jjj that the Bethe solution gives exactly the Bose- 
Einstein statistics. For c > 0, Eq.(3) allows string solu- 
tions with 



¥ - k 



+ -(/ + !- 2a)c, 



a = 1, 2, • • • , I. 



(4) 



As c — > + , the strings are squeezed to the real axis and 
hence can be treated as the super-particles. Substituting 
the squeezed-string solutions into Eq.(3) and take the 
logarithm, we readily obtain 



lk)L = 27:1) 



N 



(5) 



i=l m=l 



where 8{x) is the sign function; jj are integers (if N odd) 
or half odd integers (if N even). Notice above we have 
taken the Planck constant ft and the double mass 2M 
as our unit. Without the prefactor I, Eq.(5) is just the 
equation of the Haldane-Sharstry spin chain obtained by 
Ha and Haldane. |T^| In the continuum limit L — > oo, put 
jj -> Ij(fc). Then pt(k) + p\ (k) = dli(k)/dk as pointed 
out by Yang and Yang, jjj In this case, we derive out 
Eq.(l) from Eq.(5) in one dimension with 



flm = 2min{l,m}. 



(6) 



To verify further the validity of Eq.(l) with Eq.(6) for an 
ideal Bose gas, we study the thermodynamics of such a 
system. The density of thermal potential of this system 
can be expressed formally as 



where d is the dimension of the system and r]i(k) = 
pf(k)/pi(k) and at thermal equilibrium satisfy the re- 
lations 

ln[l + m (k)} = * [e(fc jr M] + E fh» ln t! + ^'{k)]. (10) 



After some manipulations, Eq.(10) can be converted to 
the following algebraic equations 



r,f(k) = (l + m+1 (k))(l + m _ 1 (k)), 



(11) 



with the boundary conditions lirm^oo rji(k)/l = [e(fc) 
p]/T and r] = 0. The solution of Eq.(ll) is @ 



sinh 2 [a + 1) £fa] 

mw = -^—^ i- 



sinh 2 silpL 



(12) 



From Eq.(10) we have 



^JlnCl + r?,- 1 ®) 



= \ lim [ln(l + m (k)) n'-^-^] (13) 

Z iv— >oo 1 

= -ln(l - e — 4^). 

This allows us to rewrite the density of the thermal po- 
tential as 



(14) 



It is just the expected one for a free boson gas in d dimen- 
sions. Therefore, the super-particle description Eq.(l) 
with Eq.(6) is indeed a correct one for the ideal boson 
gas in the thermodynamic limit. 

For an interacting Bose system, the total energy will be 
lifted by the repulsions among the particles. That means 
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the repulsive interaction in some sense may induce an 
exclusion effect. Suppose we introduce the interaction 
in the following way: (i)There is a one-to-one correspon- 
dence between the quasi-super-particles and the super- 
particles of the free boson system; (ii)The interactions are 
completely absorbed into the quasi-super-particles and 
induce further exclusion effect. This hypothesis defines 
the quasi-super-particle distribution as 



where 



I 



(2tt) c 



^2 fimPmik) - -fin(k). (15) 



m— 1 



Note in the above definition, we choose the exclusion 
term proportional to the occupation density n(k). 7; 
are some positive constants which measure the strength 
of the exclusion effect. A simple choice of the statistical 
factors is to take all of them to be an unique constant, 
i.e., 7; = 7. This definition is very similar to Haldane's 
but in the representation of quasi-super-particles. This 
however allows us to construct an ideal statistics without 
complete exclusion. The density of the thermodynamic 
potential in our case is still expressed as 

00 „ 

/ = £ / i\m - p] Pl {k)dk 
1=1 J 

00 

-T£ / {[ Pl (k) +p?(k)}H Pl (k) + p?(k)] (16) 
1=1 J 

-pi(k) lnp ; (fc) -tf(k) lnp?(k)}dk. 

Notice that Sp\ L (k) and Sp m (k) are not independent. 
From Eq.(15) we have 



$Pi(k) = ~Y1 + rrvy)5p m {k). 



(17) 



By minimizing the density of the thermodynamic poten- 
tial / with respect to pi(k), i.e., 5f/5pi(k) = 0, we obtain 
that 



\n[l+ m (k)} 



l[e(k)-p] 



m— 1 

With the identities 

^(/l+lm + fl-lm) = flm — Sim, I > 1 
7T flm = flm - Si 



(18) 



(19) 



Eq.(18) can still be transformed into Eq.(ll) but with 
different boundary conditions 



Vo(k) = 0, 

lim^ = ^ + 7ln^), 

I— »oo I 1 



(20) 



\nZ{k)= ^Hl + r?- 1 ^)]. 

m— 1 

The solutions of r]i(k) in the present case read 
sinh 2 (Z + l)[^%^ + ^lnZ(fc)] 



(21) 



sinh 



2T 

2\ <i(k)-n , 7 
2T ~ 2 



\nZ(k)} 



1. 



(22) 



From the / = 1 case of Eq.(18) we derive the following 
relation for the function Z(k) 



[Z{k)-l]Z-<{k)^e t ^ hl . 



(23) 



From the / — > oo case of Eq.(18) we have the relation 

oo 



1=1 



= \ Jim pn(l + ^(fc)) - n"-^— <± - N 7 lnZ(k)} (24) 

= -\n{2-Z(k)}, 

and therefore the density of the thermodynamic potential 
can be rewritten as 

f = ^kyS ln[2 ~ z[%)]dL (25) 

Since all rji{k) are non-negative, from the definition we 
know that Z{k) > 1. In addition, Z{k) < 2 to keep / to 
be real. From Eq.(23), it is easily to derive 



Z(k) - 1 



, dZ(k) = 

dp l+ 7 (l-Z- 1 (fc))' 



(26) 



With the relation J n(k)dk = —Tdf/dp we get the 
single-particle distribution function n{k) as 



n(fc) = 



Z{k) - 1 



(27r) d [2-Z(fc)][l+ 7 (l-Z-i(fc))]' 



(27) 



As 7 — 0, n(k) is just the Bose-Einstein's distribution 
function. For Z(k) — > 2, n(k) — » oo. Therefore, Eq.(15) 
can be considered as a generalization of Bose-Einstein's 
statistics. It is interesting to study the condensation be- 
havior of the present system. As in the ideal boson gas, 
condensation occurs on the k = state, if there is any. 
Note that Z(k) is a decreasing function of p. When the 
temperature is lowered, p tends to a critical value with 
which Z(0) = 2, corresponding to the singular point 
of n(k). Obviously, p(T) = 7Tln2 below the critical 
temperature T c . For a given density of particles n and 
e(k) = k 2 , T c is determined by 
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with 



Eq.(36) can be transformed into 



[z(x) - l]z\x) = 2"'e 



(29) 



Y a {k)[i + Y{k)y- 9 = 



(38) 



For x — > 0, 2 — z(x) ~ a;, the integral in Eq.(28) is di- 
vergent for <i < 2, which indicates the absence of con- 
densation at finite temperature. However, the finite- 
temperature condensation indeed occurs for d > 2. A 
simple analysis leads to dz(x)/dj > and dT c /dj < 0, 
which indicates that T c is decreasing with the increase of 
7 and T c — > for 7 — > 00. 

An interesting issue is that a proper choice of 7 m may 
even induce some exclusion statistics. For example, we 
put 7 m = mg with g a positive constant. In this case, we 
have 

oc 

Spi(k) = - £ + lmg)Sp m (k). (30) 

m— 1 

With the above relation to minimize the thermodynamic 
potential Eq.(16), we get 



ln[l +»»(*;)] = 



T 



+ J2(fim + lmg)ln[l + i 1 - 1 (k)}. 

The solutions of f]i(k) take the following form 
sinh 2 a + l)[fe + 2 1nC(fc)] 

m{k) = 1, 

where 



(31) 



sinh 



2T 



+ ?lnC(A)] 



lnC(fc) = £mln[l + ry m 1 (fc)], 



(32) 



(33) 



and the thermodynamical potential in this case reads 
1 



With the relation 



(27T)° 



In CO)*. 



(34) 



-N 



1 

~ 2 N 
e(k) - [i 



lim {ln[l +rj N (k)] 
NglnC(k)}, 



(35) 



and the distribution function n(k) reads 



n{k) 



T ainC(fc) _ 1 1 



(39) 



Eq.(38) and Eq.(39) are just the results of Haldane's ex- 
clusion statistics derived by Wu. p] 

In conclusion, by constructing the super-particle rep- 
resentation of the free boson gas, we introduce a non- 
exclusion statistics of many-particle systems, the Bosc- 
Einstein condensation of this system is also studied. Un- 
like in the free boson gas, it is found that the chem- 
ical potential below the transition temperature is lin- 
early proportional to the temperature. Haldane-Wu's 
fractional statistics is perfectly retrieved in this repre- 
sentation. This strongly suggests that with the super- 
particle representation, the non-exclusion statistics and 
the exclusion statistics can be constructed in an uniform 
way. 

This work was supported by NSFC under grant No. 
19825112. 
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